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ABSTRACT 

Satellite instruments provide global maps of surface UV irradiance by combining backscattered 
radiance measurements with radiative transfer models. The models are limited by uncertainties in input 
parameters of the atmosphere and the surface. We evaluate the effects of possible enhancements of the 
current TOMS surface UV irradiance algorithm focusing on effects of diurnal variation of cloudiness 
and improved treatment of snow/ice. The emphasis is on comparison between the results of the current 
(version 1) TOMS UV algorithm and each of the changes proposed. We evaluate different approaches 
for improved treatment of pixel average cloud attenuation, with and without snow/ice on the ground. In 
addition to treating clouds based only on the measurements at the local time of the TOMS 
observations, the results from other satellites and weather assimilation models can be used to estimate 
attenuation of the incident UV irradiance throughout the day. A new method is proposed to obtain a 
more realistic treatment of snow covered terrain. The method is based on a statistical relation between 
UV reflectivity and snow depth. The new method reduced the bias between the TOMS UV 
estimations and ground-based UV measurements for snow periods. The improved (version 2) 
algorithm will be applied to re-process the existing TOMS UV data record (since 1978) and to the 
future satellite sensors (e.g., Quik/TOMS, GOME, OMI on EOS/Aura and Triana/EPIC). 
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1. INTRODUCTION 

Global increases in UV-B fluxes from decreasing stratospheric ozone amounts 1 " 13 caused by 
anthropogenic chlorine releasing gases (mostly chlorofluorocarbons) have been an issue of public 
concern for the past 20 years because of their impact on human health, as well as terrestrial and aquatic 
ecosystems 14 20 Several satellite-based methods for estimating UV irradiance have been suggested 

4 8 21-37 

. Because of long time record and global contiguous spatial coverage, the NASA Total Ozone 
Mapping Spectrometer (TOMS) data 38 were used for estimating global trends in surface UV 
irradiance 4, 8,24 and studying global UV climatology 28,31 (especially for regions not covered by ground- 
based UV networks and over oceans 39,40 ). The TOMS UV record will continue with Ozone Monitoring 
Instrument 41 (OMI to be launched in 2004 on NASA EOS/Aura satellite) as the successor to TOMS. 

First intercomparisons with ground-based UV measurements have indicated positive bias of the current 
(version 1) TOMS UV data for many locations in NH (~10% station-average overestimation under 
snow-free conditions 31>42 " 44,76 ) and possible underestimation in the presence of snow 45 . Part of this 
bias can be attributed to the current (version 1) TOMS UV algorithm. The goal of this paper is to 
describe possible improvements to the operational TOMS UV algorithm 23, 27,31,36,41 that could help 
reducing this bias. Some of the proposed improvements will be implemented in the second version of 
the TOMS UV algorithm after extensive validation period (2002-2003). The improved (version 2) UV 
algorithm will be shared between TOMS and future GOME, TRIANA/EPIC and OMI UV products. 
The previous TOMS data record (since 1979) will be re-processed with the Version-2 UV algorithm to 
ensure homogeneity of the combined satellite global UV record for the trend analysis. 
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2. OVERVIEW OF THE CURRENT TOMS UV ALGORITHM 

The amount of ultraviolet radiation in the UVA (320nm - 400nm) and UVB (290nm - 320nm) spectral 
ranges that reach the surface of the Earth is determined by Rayleigh scattering from the molecular 
atmosphere, the absorption of ozone, scattering by clouds, and both scattering and absorption by 
aerosols 5 13 . Current TOMS UV algorithm (version 1) is based on corrections to calculated clear-sky 
UV irradiance, Eciear The calculation procedure is based on table lookup and either cloud/non- 
absorbing aerosol correction 36 or absorbing aerosol correction 27,31 (figure 1): 


F = F C 

^ cloud ^clear^-'T 


CO 


Calculation of E c i e ar in the UV range from satellite-derived spectral extraterrestrial solar irradiance 46 ’ 49 
and NASA’s Total Ozone Mapping Spectrometer (TOMS) measurements of total column ozone, 
aerosols and surface reflectivity 38 was previously described in the literature 27,31 , including estimates of 
the various error sources (Table 1). In the absence of snow, clouds, and aerosols, the effects of ozone, 
solar zenith angle, and altitude are essentially well-understood problems. Elowever, presence of 
aerosols 50 55 , clouds 5,10,13,21 ’ 36 and snow 56 ' 58 in the satellite Field of View (FOV) requires additional 
corrections. The exact correction would require complete characterization of the optical state of the 
atmosphere and the surface during the course of the day (for daily exposure calculations). The 
complete information is never available from the satellite data alone. Therefore, the correction factor 
(Ct) has to be estimated using limited information available from the single satellite measurement at 
the overpass time. Currently, TOMS UV algorithm (version 1) estimates daily exposures assuming no 
diurnal changes in Ct factor (“frozen clouds and aerosols” ). The type of correction (specific Cj 
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algorithm) is selected based on the two threshold values of the TOMS aerosol index (AI) (calculated 
from 340nm and 380nm radiances in case of Nimbus 7 TOMS and from 331nm and 360nm in case of 
Earth Probe TOMS) and Lambertian Equivalent Reflectivity (LER) (360nm or 380nm). The surface 
albedo and snow effects are estimated using the TOMS monthly Minimum Lambertian Effective 
surface Reflectivity (MLER) global databse 59,60 as described in 36,41 . 

3. IMPROVED CLOUD CORRECTION 

The main cause of daily UV irradiance variability at a given geographic location is due to clouds. Polar 
orbiting satellites allow only one daily observation at a given location. Even at the overpass time both 
cloud fraction and cloud optical thickness cannot be simultaneously derived on the basis of one TOMS 
radiance measurement. Currently, to estimate the cloud transmittance at the overpass time TOMS 
algorithm uses homogeneous cloud model embedded into Rayleigh scattering atmosphere with known 
ozone absorption and surface reflectivity, As, assuming 100% cloud cover 36,4 \ For snow-free 
conditions in the near UV spectral region As is uniformly low and can be accurately predicted from a 
climatological data base that was developed using 1 5 years of TOMS data 59 . The “effective” cloud 
optical thickness, x(to), is derived by matching the measured 380nm radiance at the overpass time, to, 
with the calculated radiance for each TOMS FOV. The same cloud model is used to calculate the 
Cj(to) as a function of x(to), As, and solar zenith angle, 9o, at all UV wavelengths, assuming spectrally 
independent x(to): 


C t OMo ) = C T (A, r(t 0 ), A s , 0 O (t 0 )) 


( 2 ) 
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The method accounts for the spectral dependence of Cj (resulted from reflection between the cloud and 
the atmosphere) as well as multiple reflections between the cloud and the surface. The second effect is 
especially important for snow covered surfaces. To calculate daily exposure, diurnal variation of 
Cr(A,t) should be estimated. For TOMS, only one observation per day is available for low- and 
midlatitude locations. Therefore, the TOMS algorithm assumes that cloud transmittance is the same as 
at overpass time: 

C T (Z,t) = C T (A,t Q ) (3) 

This may lead to erroneous UV-exposure estimates for areas with systematic changes in cloudiness 
during the day. An obvious example of such a case is for cumulus clouds, which commonly develop 
during afternoon hours (Figure 2). 

3.1 Assimilating ECMWF water content to construct a time-resolved homogeneous cloud model 

To improve daily exposure estimates additional cloudiness information per day is needed. Such 
information is available from global operational weather models and from geostationary satellites 
(GOES, Meteosat) 32 . For example, one can use vertically integrated cloud parameters (total cloud 
cover, Tc, total column water and ice content, LWC(t) provided by numerical model of the European 
Center for Medium Range Weather forecasts (ECMWF). Global ground-based, balloon borne and 
satellite weather observations are used for producing global 3-dimensional analysis every 3 hours 
(http://www.ecmwf.int/) . The spatial resolution of the Center's current model is equivalent to having 
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gridpoints separated by ~60 km around the globe. The estimated daily data volume is -20 -3 0MB /day. 
The easiest way of assimilating the model parameters into the TOMS UV processing algorithm is to 
scale TOMS effective cloud optical depth proportionally to the diurnal changes of the ECMWF model 
total column water and ice content, LWC(t): 

C T (A,t) = C T (l,T(t ),A s ,6 0 (t )) ( 4 ) 

r C) = h(‘«) Il ~ 

T e( { o) 


where Tg is ECMWF cloud optical depth, estimated from the total column water and ice content, 

assuming homogeneous cloud layer with Cl droplet size distribution. The algorithm becomes 
computationally unstable when LWC(to) (and TeOo) ) approaches zero or there are large disagreements 
between the model and TOMS estimates of the cloud optical thickness at the overpass time. In such 
cases, the following algorithm could be suggested that is computationally stable: 


T(t) = T(t 0 ) + T E (t)~T E (t 0 ) 


We note that above equation reduces to equation (4) if the TOMS and model optical thickness agree at 
the overpass time: t((q) = The scheme also works when either t((q)= 0 or TE(to)-0. If TOMS 

and ECMWF values are far apart at the overpass time (which is not uncommon), the above equation 
gives more reasonable values of optical thickness than equation (4). 
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An example of the methodology is shown in Figure 2, which shows the daily evolution of LWC(t) and 
UV radiation as measured by Brewer instrument on 10 August 1992 at Sodankyla (Finland). 
According to the model (and SYNOP observations) both Tc and LWC(t) increased during that day. 
The UV dose rates calculated with the old (Equation 3) and new (Equation 4) versions of TOMS UV 
algorithm are shown with approximately 1 0-minute time step. The dose rates calculated by the original 
algorithm (closed circles - equation 3) follow closely the clear-sky dose rates (open circles) because 
the overpass time (to= 8UTC) was almost cloud-free. As a result, the original method overestimates the 
daily UV dose by 60% as compared with the ground-based data while the new algorithm (Equation 4) 
is close to the observation. 

The proposed daily-exposure cloud algorithm was tested using combined noontime TOMS and 3-hour 
ECMWF data for summer 1992 at Sodankyla (Figure 3). The overall performance of the new 
algorithm is better than the old one: the mean difference between the TOMS UV and Brewer data 
reduced from 182 J/m 2 to 126J/m 2 during the period, at the same time correlation increased from 0.86 
to 0.89. The improvement is better than the absolute numbers indicate, since a substantial portion 
(about half) of the current bias between the TOMS and the Brewer measurements is not related to the 
TOMS cloud algorithm 31 . 

These results present a first-guess implementation using only water content time-resolved cloud data 
from the ECMWF model. Fine-tuning of the cloud algorithm will be performed to improve the 
agreement between the observed and the improved cases further. In the TOMS UV project the global 
comparability between different cloud models will be studied (i.e. ECMWF, NOAA-NCEP, NASA- 
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Goddard, NASA-ISCCP models), and improvements in the methodology will be made if needed. For 
example, we plan to compare cloud data from the ECMWF model with the similar data from 
International Satellite Cloud Climatology Project (ISCCP - http ://isccp.giss.nas a.gOv/j when both are 
available. ISCCP analysis combines satellite-measured radiances from geostationary and polar 
satellites with ice/ snow data to obtain information about clouds and the surface 61 . The pixel analysis is 
performed separately for each satellite radiance data set and the merged results are reported in the 
Stage DX data product, which has a nominal resolution of 30km and 3 hours. Using ISCCP data 
operationally for TOMS processing is currently impossible due to the large time delay in releasing 
ISCCP data products. The ISCCP data may be used for future re-processing of the global TOMS UV 
maps. We will also analyze possibility of using single satellite time-resolved cloud information in 
producing regional UV maps and overpass data sets. One possibility may be using Meteosat Second 
Generation (MSG) Climate SAF products for producing European UV maps 32 and NOAA-GOES data 
for US coverage. 


3.2 Assimilating ECMWF cloud fraction and water content with a fractional cloud model 

The overestimation of the TOMS UV data for high irradiance levels (Figure 3) may be related to the 
broken cloud effects. The ECMWF model can be combined with a fractional cloud model 62 since it 
contains the cloud cover information, T c (t). We will study a possibility to use both Tc and water 
content from the ECMWF model to improve estimation of Cj. First, the ECMWF cloud parameters 
should be adjusted using TOMS radiance measurement at the overpass time. For areas where there are 
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ground-based measurements or a region with specific biological interest, the cloud pattern data from 
high resolution satellite images (ATSR-2, AVHRR, SeaWiFs, geostationary satellites e.g., GOES) can 
be also used to estimate true cloud fraction, Tc(to), so the ECMWF model data can be adjusted. 
Second, the effective optical depth of the cloud portion of the TOMS FOV, T F (to) can be derived by 
matching the TOMS measured and pre-calculated radiances at the overpass time: 

'“('o) = TMI^ArM) + (1 - TMVSZ, (5) 

Finally, the ECMWF time resolved cloud data, normalized by the T c (t 0 ) and T F (t 0 ), can be used to 
predict diurnal variation of C T (t), similar to equation (4): 


Q (0 = 1- T c (OP - C T (X, t f (0, R s , 0 O )] 


( 6 ) 


3.3 Cloud shape effects 

The fractional cloud model described ignores the real 3 -dimensional (3D) cloud structure and some 
related cloud-radiation effects (cloud shadows, reflection from non-horizontal surfaces). We will 
estimate “cloud shape” errors and develop average Ct corrections for the TOMS FOVs using more 
realistic cloud models and Monte-Carlo radiative transfer code 63 . Figure 4 shows one possible cloud 
model, which describes fair weather cumulus cloud field. The model relates stochastic field 
characteristics with cloud amount, mean cloud diameter and aspect ratio. Based on these input 
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parameters a representation of cloud field is calculated as convolution of two-dimensional Fourier 
series with random coefficients. Calculations of the radiance at the top of the atmosphere and 
irradiance at the surface are performed using 3-D Monte-Carlo (MC) code 63 . 

Figure 5 shows the simulated normalized angular distribution of the 380nm radiance (Anisotropic 
function, AIF 64 ) backscattered from the cloud scene shown in Figure 4. The “cloud shape” error in 
surface irradiance is proportional to the AIF ratio for broken (Figure 5) and homogeneous cloud scenes 
averaged over TOMS FOV. The error is a function of assumed cloud parameters, observational 
geometry, and surface albedo. Figure 6 shows the correction factor for cloud scene shown in figure 4. 
The factor should be applied to the standard Cr value calculated using the optically equivalent (i.e. 
providing the same 380nm radiance at the satellite) homogeneous cloud model (see Equation 2). As 
expected, the factor is maximal in the solar principal plane. For this particular cloud scene and 9 0 =54° 

the factor ranges from 0.85 (<p=0°- forward reflecting) to 1.2 (<p=180° - backward reflecting, i.e. “hot 
spot”). 

Because the equatorial overpass occurs close to solar noontime, the TOMS instrument scans in a 
direction, which is approximately perpendicular to the principal plane of the sun. Figure 6 shows that 
for these directions the correction factor is much less than in the solar principal plane. However, the 
errors may be still significant for specific observational conditions (we found that the error increases 
with solar zenith angle, i.e. at high latitudes). To quantify the errors for all possible conditions the 3D 
Monte-Carlo broken cloud radiative transfer model will be combined with realistic snow BRDF 
(Bidirectional Reflection Distribution Function) to calculate the surface UV irradiance over snow with 
broken cloud conditions. Preliminary calculations have shown that broken cloud effects over snow 
differ from the uniform cloud effects over snow surface. The broken cloud model will be also 
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combined with tropospheric aerosol models to study combined cloud/ aerosol / radiation interactions. 
Before adopting either of these approaches in the version 2 TOMS UV algorithm we perform 
validation studies using radiative transfer modeling and comparisons with ground-based UV data for 
selected stations. 


4. SNOW EFFECTS 


A major problem in using satellite data to estimate UV irradiance at high latitudes arises from the 
difficulty in identifying the presence of clouds when there is snow on the ground. When TOMS views 
a scene containing ice, snow, and clouds there is no way to separate the effects of snow from clouds 
based on one reflectivity measurement. However, if the surface reflectivity (albedo), As is established 
for various conditions in a geographical region, the excess scene reflectivity can be used to estimate 
cloud transmittance, Ct, over snow surface 36,41 . Current TOMS UV algorithm uses monthly minimal 
Lambert equivalent surface reflectivity (MLER) global database 59,60 to estimate A s at 380nm. The 
algorithm also assumes that As does not change with wavelength in the UVA and UVB spectral 
regions. The assumption are in reasonable agreement with direct ground-based measurements of UV 
albedo 5,7,56 58,65 66 . MLER is a reasonable estimate of the surface albedo for either snow-free 
conditions or regions with permanent snow cover (Antarctica, Greenland). However, MLER is not a 
good estimator of actual surface albedo at mid-latitudes in spring and fall seasons when surface albedo 
varies daily depending on the presence and state of snow cover. In absence of external snow 
information the current TOMS algorithm uses a climatological snow/ice flag (probability of the 
presence of snow on a given day at a given location) to detect the presence of snow. If snow is detected 
(or likely), the algorithm first determines a snow albedo threshold (SAT). Currently the SAT is the 
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MLER value bounded by a constant value of 0.4. The value 0.4 was selected as appropriate for snow 
covered urban/suburban-populated areas containing at least moderate densities of roads, houses, and 
trees (e.g., Toronto, Moscow) 6,41 . The daily estimation of A$ is based on comparison of SAT with the 
TOMS measured reflectivity (LER) at 380nm (or 360nm). If LER is less than SAT+0.05, the cloud 
free conditions are assumed and As is set equal to LER (Cr=l). If LER is more than SAT+0.05, As is 
set equal to SAT and all additional measured reflectivity is assigned to a cloud above the snow surface. 
The algorithm proceeds to calculation of effective cloud optical thickness and Cj as described above. 
The constant threshold value 0.05 was chosen because of the difficulty in detecting thin clouds over 
snow surface. This problem exacerbates at high latitudes over areas with permanent snow cover 
(Greenland, Antarctica). For such regions, the possible error in cloud correction could exceed the error 
due to neglecting the clouds. Therefore, cloud free conditions are also assumed if the SAT value is 

more than 0.9. On average, current TOMS algorithm leads to underestimation of UV radiation in 
winter conditions 35,37,45 . 

To obtain a more realistic treatment of the albedo of snow covered terrain, a new method was 
developed . This method is based on a relation between UV reflectivity and snow depth. To establish 
such a relationship cloudless days with snow cover were selected based on the total cloud cover and 
snow depth (SD) parameters of the ECMWF ERA- 15 re-analysis dataset (covering time period 
1979-1994). These data together with the coincident and co-located Nimbus-7 TOMS reflectivity 
measurements were used to set a regression: 


Albedo = axSD u 3 


( 7 ) 
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where a is fit parameter specific to each ground pixel. As expected, this method improved the 
correspondence between the satellite-retrieved results and ground-based measurements (figure 7). 

5. AEROSOL EFFECTS 

For the purposes of estimating UV irradiance at the Earth’s surface, there are two classes of aerosols 
that must be considered: 1) aerosols that only scatter LJV radiation and 2) aerosols that both scatter and 
absorb UV radiation. The first category is included in the measured scene reflectivity (cloud optical 
thickness) and attenuates UV radiation in a manner that approximates clouds of equivalent reflectivity, 
except for the spectral dependence of the optical thickness and phase function differences. Attenuation 
values of up to 10% are common. Operational (Version 1) TOMS UV algorithm does not distinguish 
between water clouds, haze, ice clouds and non-absorbing aerosols. This could result in a few percent 
errors when a non-absorbing aerosol layer is treated as an optically thin cloud of the same reflectance. 
The error due to non-absorbing aerosols could be both positive and negative depending on the 
observational geometry (phase function effect) and could be treated as a random error, rather than the 

constant bias. 

On the other hand, a significantly larger overestimation occurs if absorbing aerosol is treated as non- 
absorbing or clouds 27 . Radiative transfer calculations with different absorbing aerosols models 
(industrial, smoke, dust) have shown a positive bias of the TOMS UV data if absorbing aerosols are 
treated as non-absorbing clouds. The bias is always positive (TOMS overestimates surface UV 
irradiance) and is proportional to the aerosol absorption optical thickness, which is a product of 
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extinction optical thickness and single scattering co-albedo (ratio of aerosol column absorption to the 
sum of absorption and scattering). The properties of absorbing aerosols vary depending on location, 
season, pollution and aerosol transport. Therefore, to characterize the bias for particular region a joint 
statistics of aerosol optical thickness and UV single scattering albedo is required. Not until recently 
such joint statistics becomes available from the inversions of the AERONET sun and sky photometer 
measurements 70 ’ 75 . The AERONET retrievals of single scattering albedo are currently restricted to the 
visible spectral region (shortest wavelength 440nm) and the retrieval uncertainties are large for a 
relatively clean locations (aerosol optical thickness less than 0.5) covering most of the globe. 
Therefore, the problem of extrapolating AERONET data from the visible to the UV spectral region 
remains a big challenge. 

Under certain favorable conditions absorbing aerosol plumes transported in a free troposphere (dust 
and biomass burning smoke plumes at high altitudes) can be detected directly using the TOMS Aerosol 
Index (AI) data 68 ’ 69 . For such cases (mostly in tropical regions) the cloud correction is replaced by 
absorbing aerosol correction (AAC) algorithm 27,31 . The AAC algorithm accounts for larger attenuation 
of UV irradiance by absorbing aerosols compare to clouds/non-absorbing aerosols of the same 
reflectivity (see figure 1). Currently, the AI data are used directly to correct Ej ear ( 1 ) for attenuation by 
absorbing aerosols: 


^ aerosol _ g (H A ) AI 

F 

^ Clear 


( 8 ) 
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where conversion factor g is a function of aerosol height, Ha and aerosol type. For the AAC method 
the major problem arises from uncertainty in aerosol plume height. Current TOMS AAC algorithm 
assumes the nominal height of 3km for plumes of desert dust and biomass burning smoke in the tropics 
27,31 . The uncertainty in the actual aerosol height is included in the error budget of the TOMS UV 
products as shown in table 2 31 . In the second version of the TOMS UV algorithm the Ha could be 
estimated using the GSFC data assimilation model or other sources. 

The AAC technique does not have the sensitivity to detect weakly absorbing aerosols in the boundary 
layer. This is the main reason of overestimating surface UV under cloud-free conditions for industrial 
regions. The treatment of aerosol attenuation of UV irradiance might have been improved by using 
newly developed TOMS aerosol products: optical depth and single-scattering albedo for dust, smoke, 
and sulfate aerosols 69 . However, the method requires a-priori knowledge of the aerosol vertical 
distribution and information about aerosol parameters (either size distribution or refractive index). 
Complete aerosol retrieval is possible by combining TOMS measurements with ground-based active 
(lidar) or passive (sun and sky photometer) measurements. Combining the TOMS data with sun- 
photometer (AERONET network) and lidar data appears to be the most perspective way of improving 
aerosol correction algorithm for regions affected by smoke from biomass burning or desert dust. 
However, characterization of the aerosol properties (size distribution and single-scattering albedo) in 
an atmospheric column from the sky photometer measurements becomes increasingly challenging for 
small amounts of absorbing aerosols. In addition, presence of broken clouds adds to the uncertainty of 
both ground-based and TOMS aerosol retrievals. Thus, improving characterization of the average 
transmittance through the inhomogeneous scenes including mixtures of broken clouds with absorbing 
aerosols remains the most challenging problem in satellite UV estimation. 
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6. CONCLUSIONS 

We evaluated the effects of possible enhancements of the current (version 1) TOMS surface UV 
irradiance algorithm. The enhancements include more detailed treatment of tropospheric aerosols, 
effects of diurnal variation of cloudiness and improved treatment of snow/ice. Some of the proposed 
improvements will be implemented in the second version of the TOMS UV algorithm after extensive 
validation period (2002-2003). Validation of spatial average UV irradiance (satellite) with temporal 
average UV data (ground station) under broken cloud conditions would require implementation of 
special sub-satellite UV validation campaigns and optimal ground UV validation strategy. 
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Table 1 Expected errors in estimated surface flux at 310nm 


IdUI 

Error source 

SZA=0° 

SZA=50° 

sza=70° 

sza=86° 

Total Ozone (2% rms) 

1.3% 

2% 

3% 

< 3% (Umkehr) 

Surface reflectivity 
snow/ice- free, 1% rms 
(* snow/ice, 5% rms) 

0.4% (3%) 

0.4% (3%) 

0.4% (3%) 

0.5% (3%) 

Strat profile (1/4 th of high- 
low difference) 

0.5% 

0.7% 

0.3% 

>30% (Umkehr effect) 

Trop profile (5 DU rms) 

0.5% 

0.5% 

0.2% 

0.2% 

Total rms 

2% (3.6%) 

2.8% (4%) 

3.7%( 4.7%) 

>30% 


The table concerns only cloud- and aerosol-free cases and gives lower limits of the uncertainties for horizontally homogeneous 
scenes over OMI footprint; (2) Numbers in parenthesis apply to snow/ice conditions for horizontally homogeneous surfaces wit 
high surface albedo, R s > 0.7 (Antarctic, Greenland). Over continents, a more realistic rms for snow/ice al * edo ™ 0 ^*** 1 ^ 
10-15%; (3) Errors do not include a ± 3 % uncertainty in absolute flux attributable to extra-terrestrial solar flux, (} ty 

in long-term trend is 1%/decade (Herman et al. 1996), caused primarily by total ozone uncertainties; (5)The flux at 3 lOnm has 
roughly the same ozone response as the erythema-weighted flux. Uncertainties in spectral flux at other wave engt s for 0 
can be approximately scaled relative to 310nm by using the following factors: 2 (at 305nm), 4 (at 300nm), and 6 (at 295nm). 

6 0 > 50* the uncertainties in spectral flux can not be scaled because of the Utnkehr effect. 
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Table 2 Expected Errors in UV-Irradiance Attenuation (AAC) algorithm due to uncertainty in 
estimated aerosol height, AH, for smoke and dust 31 . The error is proportional to both AH and aerosol 
absorption (single scattering co-albedo, 1-co). The errors are shown in Table 2 for AH=0.5km. The 
largest error corresponds to the large-particle dust, which is also strongly absorbing . 


Model/Parameter 

CO 

k 

S (km' 1 ) 

E f (%) 

Smoke Cl 

0.92 

0.25 

12 

1.5 

Smoke C2 

0.84 

0.37 

40 

7.5 

Dust D 1 

0.90 

0.28 

12 

1.7 

Dust D2 

0.72 

0.50 

43 

11 

Dust D3 

0.63 

0.57 

55 

16 




29 


Figure captions 

Figure 1. TOMS version 1 UV algorithm processes diagram 

Figure 2. Comparison of original and improved cloudiness treatment in the UV algorithm. Liquid and 
ice water data of ECMWF (solid line) is used (assuming total overcast case) in the improved case. The 
improved method (Eq.4) shows better agreement with the Brewer daily observations (diamonds), 
whereas the original TOMS algorithm (Eq.3) gives an overestimate of the daily dose -60% for this 
particular case. 

Figure 3. Comparison of daily Erythemal (CIE) UV doses [J/m 2 ] calculated with original (Eq.3) and 
improved (Eq.4) cloudiness treatment in the TOMS UV algorithm with ground based Brewer 
observations at Sodankyla from 8 July to 1 September, 1992. 

Figure 4. Fragment of he broken cloud field as model input. The model is based on the normal 
random (Gaussian) field with a fixed lower boundary [Geogdzhaev et al., 1997], Spatially 
inhomogeneous cloud structure is described by 3D array of cells each with prescribed cloud properties. 
The dimensions of each cell should be sufficiently small compare to the photons free path. In most 
cases 50 to 1 00 meters cells were found to be sufficient in modeling TOMS and OMI FOVs. The 
dimensions of the calculated field are similar to a single OMI FOV (10km by 20km). Cloud cover 0.5, 
aspect ratio 1, scattering coefficient 50 km' 1 and cloud average diameter 1km. 

Figure 5. Anisotropic function AIF(8 0 ,G,<p) of the broken cloud scene (figure 4) for0 o =54 0 as a 

function of the satellite vertical angle, 0 (giving by distance from the center of the figure - nadir 
direction) and solar azimuthal angle, <p (given by polar angle: forward reflecting in on the right and 
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backward reflecting — on the left). The AIF is defined as the ratio of the equivalent Lambertian flux to 
the actual reflected flux: AIF = ttL(9 0 ,9,<p)/ M(9 0 ) [Suttles et al. 1988], AIF <1 (shown by green color) 
means that the measured radiance is less than Lambertian and AIF>1 (red colors) means that the 
radiance for broken cloud scene is greater than Lambertian. radiance. 

Figure 6. Example of a correction factor for broken cloud scene (figure 4), which should be applied to 
the satellite UV data calculated with TOMS method (i.e. using homogeneous plane-parallel cloud 
model). Depending on the satellite viewing direction (explained in figure 5) the correction factor 
ranges from 0.85 to 1 .2. Cloud anisotropy is much less in the plane perpendicular to the solar principal 
plane than in the solar principal plane. Surface reflectivity 5%. Solar zenith angle 54°. The equivalent 
optical thickness of the homogeneous plane parallel cloud layer is close to 5. 

Figure 7. The effect of the new snow albedo treatment on the computed surface UV. “TOMS original 
UV” is based on MLER 59,60 whereas “TOMS modified UV” is based on the snow albedo regression 
with ERA-1 5 snow depth data, and coincident and co-located Nimbus-7 TOMS reflectivity 
measurements 67 . Top: Comparison at Sodankyla, Finland; bottom: Churchill, Canada 
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Figure 1 TOMS version 1 UV algorithm processes diagram 







32 



Figure 2. Comparison of original and improved cloudiness treatment in the UV algorithm. Liquid and 
ice water data of ECMWF (solid line) is used (assuming total overcast case) in the improved case. The 
improved method (Eq.4) shows better agreement with the Brewer daily observations (diamonds), 
whereas the original TOMS algorithm (Eq.3) gives an overestimate of the daily dose -60% for this 


particular case. 
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TOMS original UV TOMS modified UV 



Figure 3. Comparison of daily Erythemal (CIE) UV doses [J/m ] calculated with original (Eq.3) and 
improved (Eq.4) cloudiness treatment in the TOMS UV algorithm with ground based Brewer 
observations at Sodankyla from 8 July to 1 September, 1992. 
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Figure 5. Anisotropic function AIF(0 o , 0, q>) of the broken cloud scene (figure 4) for9o-54 0 as a 

function of the satellite vertical angle, 9 (giving by distance from the center of the figure - nadir 
direction) and solar azimuthal angle, (p (given by polar angle: forward reflecting in on the right and 
backward reflecting — on the left). The AIF is defined as the ratio of the equivalent Lambertian flux to 
the actual reflected flux: AIF = nL{6 a ,9,(p)l M(0 O ) [Suttles et al. 1988]. AIF <1 (shown by green color) 
means that the measured radiance is less than Lambertian and AIF>1 (red colors) means that the 


radiance for broken cloud scene is greater than Lambertian.radiance. 
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Figure 6. Example of a correction factor for broken cloud scene (figure 4), which should be applied to 
the satellite UV data calculated with TOMS method (i.e. using homogeneous plane-parallel cloud 
model). Depending on the satellite viewing direction (explained in figure 5) the correction factor 
ranges from 0.85 to 1.2. Cloud anisotropy is much less in the plane perpendicular to the solar principal 
plane than in the solar principal plane. Surface reflectivity 5%. Solar zenith angle 54°. The equivalent 
optical thickness of the homogeneous plane parallel cloud layer is close to 5. 
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Figure 7. The effect of the new snow albedo treatment on the computed surface UV. “TOMS original 
UV” is based on MLER 59,60 whereas “TOMS modified UV” is based on the snow albedo regression 
with ERA- 15 snow depth data, and coincident and co-located Nimbus-7 TOMS reflectivity 
measurements 67 . Top: Comparison at Sodankyia, Finland; bottom: Churchill, Canada 






Popular Summary: 

Version 2 TOMS UV algorithm: problems and enhancements 

Nickolay Krotkov 1,2 , Jay Herman 2 , P.K.Bhartia 2 , Colin Sefitor 3 , Antti Arola 4 , Jussi 
Kaurola 4 , Lasse Koskinen 4 , S. Kalliskota 4 , Petteri Taalas 4 , 1. Geogdzhaev 5 


Satellite instruments provide global maps of surface UV irradiance by combining 
backscattered radiance measurements with radiative transfer models. The models are 
limited by uncertainties in input parameters of the atmosphere and the surface. We 
evaluate the effects of possible enhancements of the current TOMS surface UV irradiance 
algorithm focusing on effects of diurnal variation of cloudiness and improved treatment 
of snow/ice. The emphasis is on comparison between the results of the current (version 1) 
TOMS UV algorithm and each of the changes proposed. We evaluate different 
approaches for improved treatment of pixel average cloud attenuation, with and without 
snow/ice on the ground. In addition to treating clouds based only on the measurements at 
the local time of the TOMS observations, the results from other satellites and weather 
assimilation models can be used to estimate attenuation of the incident UV irradiance 
throughout the day. A new method is proposed to obtain a more realistic treatment of 
snow covered terrain. The method is based on a statistical relation between UV 
reflectivity and snow depth. The new method reduced the bias between the TOMS UV 
estimations and ground-based UV measurements for snow periods. The improved 
(version 2) algorithm will be applied to re-process the existing TOMS UV data record 
(since 1978) and to the future satellite sensors (e.g., Quik/TOMS, GOME, OMI on 


EOS/Aura and Triana/EPIC). 



